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Abstract 


The conversion of particle- energy spectra into a linear energy 
transfer (LET) distribution is a guide in assessing biologically sig- 
nificant components. The mapping of LET to energy is triple valued 
and can be defined only on open subintervals. A well-defined numeri- 
cal procedure is found to allow generation of LET spectra on the open 
subintervals that are integrable in spite of their singular nature. 


Introduction 

A customary assumption is that the biological re- 
sponse (or risk) to radiation exposure is related to the 
energy absorbed per unit mass within a macroscopic 
volume of the biological material. This implies a cer- 
tain uniformity in the energy deposited (which is true 
for the radiation of low linear energy transfer (LET)) 
such as from X-ray, 7 -ray, and /Tray exposures. As 
experience was gained from exposures with high-LET 
radiation, the nonuniformity of the energy deposit 
was clearly an important determinant of biological 
response, and LET-dependent relative biological ef- 
fectiveness factors came under intense study and pro- 
vided the basis for the LET-dependent quality factor 
for use in radiation protection (ref. 1). The micro- 
dosimeter concept was developed as an instrument to 
measure such LET fluctuations. 

Although the concept of an LET-dependent rel- 
ative biological effectiveness (RBE) was useful for 
charged particles of relatively low kinetic energy, con- 
cerns for the biological effects of the high charge and 
energy ions (HZE) in deep space result from the lat- 
eral spread of the energy deposit that is larger than 
a typical cell nucleus (ref. 2). In this case, LET is 
not sufficient to define an RBE, but the track width 
also plays an important role. This was the origin of 
the limited LET concepts (refs. 1 and 3) as related 
to the lateral spread of the energy deposit. The im- 
portance of such track-structure effects on biologi- 
cal response has been demonstrated for a number of 
biological systems (ref. 4). 

Even though LET as a biological-response indica- 
tor is most applicable to low charge Z and low-energy 
particles, LET is considered as a rough indicator for 
HZE particles as well. For this reason, we are in- 
terested in the generation of differential LET distri- 
butions as a guide to identify biologically significant 
components. 

Differential LET Spectra 

In radiobiology, relat ing biological response to the 
LET of the radiation environment is a tradition. For 
example, RBE and the quality factor Q are generally 


taken as being related to LET. As a consequence, 
the concept of LET spectra has played a role in 
estimating biological response. Unfortunately, this 
concept is most useful if the flux ( <f>i{L ) dL) of 
particles with LET (that is, L) between L and L + dL 
is known. This is generally found by knowing the 
energy flux {4>e{ e ) dE) of particles with energy (E) 
between E and E + dE , where L is known as a 
function of E so that 


0L (L) = 


dL 

dE 


4>e(E) 


(l) 


where the prefactor on the right-hand side is the Ja- 
cobian between the E and L spaces. The difficulty 
with this approach is that dL/dE = 0 at the max- 
ima and minima of the LET curve; 4>p{L) must be 
replaced by the sum over the various branch func- 
tions as 
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where Eq is the energy of each branch associated 
with L. That is, for all values of E B) 

L = L (E B ) (3) 


Clearly, <Pl(E) does not exist for every value of L but 
is defined on open intervals not containing values for 
which dL/dE = 0. Furthermore, <}>l(L) is unbound 
on the open subintervals over which it is defined, even 
though 4>l{L ) is integrable over its domain. From 
the above arguments, enough challenges obviously 
exist in finding a representation for </>/,(!/), especially 
as some data set to be used by others in specific 
applications. This problem is simplified since L(E) 
has but one maximum and one minimum other than 
at zero energy. Furthermore, L(E) has continuous 
second derivatives allowing a Taylor series expansion 
as 


L{E)*L{E B )+ l -L"{E B ){E-E B ) 2 (4) 



in the neighborhood of the branch limits. However, 
one finds 
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in the neighborhood of the branch points, where the 
subscript B denotes evaluation at the branch limit. 

We implement the above considerations as fol- 
lows. The LET is defined over a numerical grid given 
by the sequence {E{\. The maximum and mini- 
mum branch points are found at dL/dE = 0 and 
are noted by E mQX and E in j n , respectively. The se- 
quence {E{)i is defined as those values oT E $ less 
than EVnax , with the main sequence {E t } m being de- 
fined by Emax < E\ < E m \ n and the sequence {Ej}jj 
defined by E m \ n < E^. The three branch functions 
are then represented by 


{^L^B 
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where B denotes one of the three branches (that 
is, B = L, m, or H). Giving a table of values 
Ei )} b in order to reconstruct an adequate 
representation of the function over each branch is not 
sufficient because 4>U is unbound near the. branch 
limits and an extrapolation into the neighborhood 
of the branch limit must be provided. If {E^} is 
sufficiently close, then the LET spectrum may be 
represented as 


^L(L)^<PLi(\Li-L D \/\L-L B \) V2 (7) 


where is the nearest grid value to the branch 

limit L& in the appropriate domain. Thus, the 
data set required to reconstruct the LET spectrum 
is the branch limit values of E max , E m j n , L max , and 
E m \w and the sequences {i?*}, {Li}, and {^Ei}- The 
numerical values of the above parameters depend on 
the charge and mass of the particles of the field. 
Thus, E t nax, E m \ n , L max , and L mm must be specified 
for each ion type in the radiation field. 

Stopping Powers 

The energy imparted to the medium by a passing 
charged particle is related to the stopping powers 
(the energy loss per unit distance traveled) and is 
known as the linear energy transfer (LET), where 
L = S(E). To calculate the spectrum of energy 
deposit, we must evaluate the LET spectrum given 
by equation (1). Clearly, a method of generating 
dS/dE as a function of E is required. As explained 
elsewhere (ref. 5), the stopping powers used are 


taken from the fitted curves for protons and alpha 
particles of Ziegler at low energy (ref. 6), and they 
are extrapolated to the Bethe formula at high energy. 
The stopping power of ions with a charge greater 
than 2 is scaled using the effective charge formalism 
and the alpha particle stopping power. The stopping 
power is stored in a data array for 15 ion charge 
values (1 to 92) over an energy grid of 60 points 
between 10 keV and 50 GeV. Extrapolation to lower 
energy and higher energy is accomplished by making 
the stopping power and its first derivative continuous 
at the boundaries of the energy grid. The numerical 

n E, Z ) j to minimize 

numerical errors. 

The scaled array is fitted by a two-dimensional 
spline with appropriate boundary conditions from 
which the derivative of S(E) with respect to E is 
found. The derivatives, which are shown in figure 1, 
vanish at the LET maximum at low energy and the 
LET minimum at high energy because of relativistic 
effects. The energies associated with the maximum 
are shown in figure 2 along with the corresponding 
energies at the minimum. Figure 3 shows the cor- 
responding values of L max and L m \ n . These provide 
the basic data for constructing differential LET (that 
is, L ) spectra. 


values are scaled by In 5(1 



F, MeV 


Figure 1. Derivatives of LET for several ions. 


Equilibrium Spectra 

Regardless of whether the source of ions is at the 
boundary or an interior volume source, the fluence 
approaches a characteristic equilibrium spectrum 
that depends on the composition of the medium. The 
equilibrium spectrum for volume sources has been 
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Figure 2. Branch limits as a function of atomic number. 



Figure 3. Values of LET at branch limits. 


given elsewhere (ref. 7). The fluence within a medium 
at a distance x from the boundary is 

cf>(x,E) = ?^<f>(0,E x ) (8) 

where 0(0, E) is the fluence at the boundary, 
E x — R~^[R(E) + x] is the energy at the boundary, 
and R(E) is the range of the particle. The <p(x, E) 
term of equation (8) results from (ref. 7) 


which is the equilibrium spectral limit. The LET 
spectra for x - 0, 0.5, 1, 3, and 5 cm of tissue are 
shown in figure 4 for three large solar events of solar 
cycles 19 and 20. The nuclear reaction products 
are neglected (refs. 7 and 8). One may observe in 
figure 4 that the equilibrium spectrum is achieved 
rather quickly above 40 MeV/cm corresponding to 
12 MeV protons. Equilibrium is achieved above 
10 MeV/cm for the 5-cm depth corresponding to 60 
to 70 MeV protons. The limit of a pure equilibrium 
spectrum (for particle type j ) at all energies is given 

= (U) 

where c is a constant. The equilibrium differential 
LET spectra for c = 1 are shown in figure 5. 

Integral LET Spectra 

The integral LET spectrum is given as 

/*Lrnax i . 

$( >L ) = J l 4>l (E ) dL (12) 

which may be related to integral energy spectra as 

$(> L) = $(> Ei) - $(> E 2 ) + $(> £3) (13) 

where E\, £2, and £3 are the three roots (branch 
functions) of 

S{E) = L (14) 

The three branch functions of equation (13) are 
shown in figure 6 for a range of L values. The 
integral LET spectra for the three solar flares shown 
previously are shown in figure 7. To gain perspective, 
we also evaluate the integral equilibrium spectrum of 
equation (11) as 

<S>(>E) = J e = R(0 - R (E) (15) 

where we arbitrarily take the variable £ to be 10 GeV. 
The integral LET spectrum for £ > £3 is then 

*(>L) = R (£2) - R (Ei) + R (0 - R (E 3 ) (16) 


dx 



<j> (x, E) = Nuclear gains and losses 


( 9 ) 


where the nuclear effects are assumed to be 0. We 
note that when x » £(£), then E x « R X (x) and 


sf/r 1 ^)] 4> [0, R~ l (x) 

S(E) 


(10) 


If we examine only the high-LET region, then 

*{> L) = R{E2)-R{Ei) (17) 

Results from equation (17) are shown in figure 8. 
Thus, the integral LET spectra are characterized by 
the main branch of the LET curve as 

$(>£)«£ (£ 2 ) 


<f> (x, £) « 


(18) 
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Figure 4. Differential LET spectra for three solar flares of 
cycles 19 and 20. 



Figure 5. Equilibrium differential LET spectra for several 
ions. 



Figure 6. The three branch functions for protons. 


except at the highest LET values where 

^(> L)^(E 2 -E 1 )/S(E m ax) (19) 

The integral LET spectra associated with the three 
solar events in figure 4 are shown in figure 7 and 
should be compared in shape to the equilibrium curve 
for protons in figure 8. 

As an application of the present procedure* the 
differential LET spectra for the solar minimum en- 
vironment are shown in figure 9 at three depths 
in tissue-equivalent material behind a slab of alu- 
minum shielding of 5 g/cm 2 . The reduction of high- 
LET components at larger tissue depths increases the 
low-LET component distributions as shown. The 
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Figure 7. Integral LET spectra for three solar events of 
cycles 19 and 20. 



Figure 8. Equilibrium integral LET spectra for several ions. 



Figure 9. Differential LET spectra at solar minimum in tissue 
shielded by 5 g/cm 2 of aluminum. 



Figure 10. Integral LET spectra at solar minimum in tissue 
shielded by 5 g/cm 2 of aluminum. 
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corresponding integral LET spectra are shown in 
figure 10. 

Concluding Remarks 

The problem of generating linear energy transfer 
(LET) spectra has been clearly delineated. Although 
integral LET spectra are continuous functions, the 
differential LET spectra are defined only over open 
subintervals and are unbound near the open end 
points. An accurate method of constructing the 
spectra over their domain has been given. These 
methods will be useful in analyzing amplications of 
shielding for protection of biological and electronic 
equipment. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
October 19, 1992 
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